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INTRODUCTION 



A poor thermal conductivity of all high-temperature superconductors makes them prone to local overheating when- 
ever high current density is applied to them. Several works used short-pulse techniques to minimize effect of heating 
in stacks of intrinsic tunneling junctions (mesas) (ITJ) 0, 0, 0, 0, IE • There were few attempts to measure directly 
^ ^ ^. the temperature of the stack utilizing the temperature-dependent c-axis resistance 0, Q- In such experiments, an 
^ ■ extra stack placed on top of- or next to |^ the major one served as thermometer. Such a thermometer could not 
c/^ \ yet correctly measure temperature under the bias lead situated at relatively large distance from the thermometer. 

■ Due to larger current density in the area of bias contact, especially if there were many "extra" junctions underneath 
2 \ it because of technical problems 0, such contacts is the major source of heating and the temperature there may be 
C . even higher than deduced in Refs. j3, 

On the other hand, estimations of overheating, based on a simple theory and numerical simulations showed that 
^ . the self-heating is proportional to a characteristic in-plane size of the stack 9j. Experimental data for small Bi2212 
Q ■ mesas were in qualitative agreement with the calculations. Estimations of self- heating in Bi2212 mesas with different 
O \ sizes demonstrated that self- heating can effectively be obviated in small mesa structures 9]. 

To better address the heating issue, we directly measure temperature of the stack by using a micron-sized thermo- 
couple which is in direct thermal- and electrical contact to the stack, exactly in the place where the bias current is 
injected. We believe that the temperature measured in this geometry is closest to the real one inside the stack. 
^ ■ 

^ . ii. samples and measurements 

On 
O 

^Nf^ ■ We have performed simultaneous temperature- and I-V measurements on a number of Bi-based single crystals, both 
\ 2201 and 2212, with three different sizes of stacks on every crystal, 4x4, 6x6, and 8x8 fim^ in area, and containing from 

■ 15 to 30 ITJ, i.e. about 200-400 A high. Those were wet-etched (in saturated solution of ethylenediaminetetraacetic 
C3 ■ acid (EDTA)) [l^. 
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FIG. 1: Schematic layout of simultaneous temperature- and I-V- measurements. Silver thin- film contacts are used for current 
(I) and voltage (V), while the bismuth one is used to monitor the temperature of the stack (red) by measuring voltage V2. 
Black arrows represent current flow. 
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Fig. 1 schematically shows a layout of contacts in the vicinity of stack (red). All the electrical contacts were 
patterned from metal thin films using photolithography and wet- and dry etching. The current ("I") and potential 
("V") leads (light-blue) are formed by just one continuous silver thin film (200 nm thick). A 200-nm-thick bismuth 
thin film (magenta) is making a counter- electrode of the Ag-Bi thermocouple. Simultaneously with the ones deposited 
onto the BSCCO-stacks, the silver and bismuth thin films were deposited on a separate glass substrate [Ojf. The latter 
formed an equivalent thermocouple which was calibrated using small diode thermometers, in a separate experiment. 
Rather large thermo-e.m.f. of such a thermocouple, about 50 /iV/K at 100 — 300 K (see Fig. 2), assures high accuracy 
of temperature measurements. 



Due to much larger voltages V ~ 1 V appearing on the potential Ag-contact, Bi-contact with V2 < 10 meV is 
important to accurately align relative to the stack and other contacts. Misalignment results in picking some voltage 
drop (very small compared with V, but essential compared to V2) along the current lead leading to highly asymmetric 
y2(/)-dependence (Fig. 3a). When the current and potential contacts were swapped, V2{I or V) became much more 
symmetric. Fig. 3b, implying that the Bi-contact was originally slightly shifted towards I-contact. This allowed us 
to perform additional check to make sure that we measure T exactly on the top of the stack. From the symmetry 
of V2{I) curves and their independence on configuration of voltage and current leads for a given stack, it is possible 
to say whether the Bi-contact was placed correctly. Measurements of all voltages and currents were done with help 
of instrumentation (pre)-amplifiers AMPOl (Fabr Analog Devices, input noise 5 uV/a/Hz) and 16-bits-20 kHz data 
acquisition card (PCI-16MI0 from National Instruments) equipped with sample- and- hold input stage for simultaneous 
sampling of data. Normally, a slow bias sweep (0.05 Hz) was used to measure both the temperature and I-V's and 
somewhat faster (0.5 Hz) sweep for just I-V's, at the same ambient temperature. Totally, 2000 — 3000 data points 
are normally logged in the former-, and 10000 in the latter case. The calibration curve shown in Fig. 2 was used to 
calculate T at each I — V point. 



III. RESULTS 



In all cases, the temperature of the stack Tg can easily rise as high as 200 — 300 K above the ambient temperature 
To, see Figs. 4-6. This means that the humps which are seen in intrinsic tunneling spectra at high bias, are much 
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FIG. 3: V{I) and V2(I) for two configurations of potential and current leads in the case when the Bi-contact was not accurately 
aligned, a) V- and I- contacts as scetched in Fig. 1; b) - with V- and I- contacts swapped. 
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FIG. 4: dI/dV{V) and temperature of the Bi2201-Lao.4 (optimal doping) stacks. Blue empty dots correspond to Ts = Tc. 
Left: 6x6-, right: 4x4 /im^ in area; 



likely due to overheating, as correctly stated in the Comment by V.N. Zavaritsky [TlL IT^. 

Fig. 4 shows dI/dV{V) and temperature of the Bi2201-Lao.4 (optimal doping) stacks [0]J. The former curves 
show all characteristic features of a typical intrinsic-tunneling spectrum. There are superconducting-gap (SO) peaks, 
then, at higher voltages - " d ips" A3i , eventually followed by " humps" . The latter two were usually attributed to the 
pseudogap (PG) 111 111 111 El 



Our temperature measurements have shown that the SG-peaks develop at voltages, corresponding, to higher than 
ambient temperature (yet at T < Tc), that agrees with what was qualitatively shown earlier 0, H 0, S US? S 0- 
However, the PG humps always lie at T > Tc, see Figs. 4 and 5. Big empty marks indicate voltages at which Tg ^ Tc. 
In Fig. 5, the vertical black lines indicate positions of the SG peaks which are not well resolved in the plot. For "8x8" 
stack the lines correspond to ~ 18 — 19 K, for "6x6" -~17 — 18K, and for "4x4" - ^ 14 K, confirming that for small 
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FIG. 5: dl /dV{V) and Ts{V) for three stacks of different areas, 8x8 (red), 6x6 (green), and 4x4 jiw? (blue) prepared on one 
and the same single crystal of Bi2201-Lao.2 (overdoped). All stacks had 12 IJJ. Triangular marks indicate the voltages at which 
Ts Tc ~ 24 K. Vertical black lines indicate the positions of superconducting peaks (although not well seen for this sample). 
For "8x8" stack, they correspond to ^ 18 - 19 K, for "6x6" - ^ 17 - 18 K, and for "4x4" - 14 K, confirming that for small 
sizes the overheating is smaller |^]. 

sizes the overheating is smaller 9j. The SG is accordingly smaller at higher temperature that explains the difference 
in SG's between samples with the same number of IJJ. 
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FIG. 6: dI/dV{y) and temperature of the HgBr2-Bi2212 8x8 /xm^ stack. Empty triangles correspond to Ts = Tc ^ 63 K. Note 
flat parts of Ts(y)-curves at —1 < V < 1 Volts. Within this range, the overheating is less than 5-6 K, see Fig. 7. 



To reduce internal heating in stacks of IJJ, one can use the so called intercalated samples, where relatively big 
chemically neutral HgBr2-molecules are interspersed inside Bi2212 single crystals 18]. The intercalated molecules 
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were found to locate themselves in between adjacent BiO layers, thus pulling apart superconducting Cu02-layers. 
This results in dramatic reduction of the c-axis critical current density (typically 10-30 times) and increase of the 
tunneling resistance, and in very much anticipated decrease of the internal heating These intercalated samples 
allowed us to reach the normal-state tunneling parts without visible signatures of internal heating (back-bending). 
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FIG. 7: /(y)-curves and temperature of the HgBr2-Bi2212 8x8 /xm^ stack, the same as in Fig. 6. The Ts{V) multi-branch 
curves were calculated from the (experimentally measured) average Ts((3)-dependence, where Q is the dissipated power. Note 
that, despite a large number of IJJ in this particular stack, the overheating does not exceed 5 — 6 K within 1 V of bias voltage. 



However, from Fig. 6 we see that internal heating is still a problem at relatively high bias even in those samples. 
Ts ^ Tc at voltages which correspond to the positions of the dips. 




V(mV) 

FIG. 8: Multi-branch hysteretic /(l/)-curves and temperature of the Bi2212 6x6 fim^ stack. For this particular stack, the 
overheating does not exceed 20 K within ±1 V of bias voltage. 
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Overheating at biases corresponding to the hysteretic multi-branch I{V) characteristics is nevertheless relatively 
small. For the intercalated sample, < 5 — 6 K at most, despite the number of junctions is big (33, compared to 
a "normal" for our experiments number of 10-15). For the pristine Bi2212- single crystals, this overheating is about 
4-5 times larger, see Fig. 8. 

It is important to take into account much likely different thermal conductivities of pristine and the intercalated 
samples, while comparing the results for them. It was noted, that the heating rate, dT/dQ^ is different, dT/dQ ~ 
30 K/mW for the pristine, and ~70 K/mW for the intercalated sample. dT/dQ was found to be almost constant for 
high temperatures and has tendency to somewhat increase on cooling. 



10 



tr48#1(6)_20K=const_d3K 



G 6 
> 

e 4 




12 



10 



tr48#*(8L1 2K=const_pm5K 



G 

73 



1 1 1 1 1 1 1 

\ Bi2201-La^j^ 


\k 8x8 ^im^ /f 




200 K 






80 K \\ V 






/ 140 K 


40 Ky 




1 1 1 


V^12K 

1 1 1 1 



-800 -400 400 

V(mV) 



800 -800 -400 400 

V(mV) 



800 



FIG. 9: Bi2201-Lax. The "heating- free" tunneling-conductance curves deduced from the "pool" of numerous "3D" I — V — Ts 
data points taken at about 30 — 40 bath temperatures in the range 4 — 220 K. The (t{V) — -points in these plots 

correspond to temperatures only different to within ±1.5 K from the temperatures shown along with the curves (with one 
exception for T^o = 12 K to the right, for which = ±2.5 K). Note also that (riV) oc \V\^ ^ with a <1 for low temperatures. 
Lines are guides to the eye. 



The detailed measurements of the temperature of the stacks at different ambient temperatures Tq allow us to 
deduce the unique, " heating- free" liV)-^ or dl / dV {V)-(ZMiYe^. First, we put the measurements for a given sample 
{I — V — Tg-triples) at all Tq, into a single file. Then, we pick out those I — V pairs which correspond to a specified 
Tgo and JT, i.e. "allowed" range of overheating. In the following, we choose 5T = ±1.5 K, counted from the chosen 
Tso- Figs. 9 and 10 show results of such a selection. 

Figs. 9 and 10 show that, after the heating is properly taken into account, the resulting tunneling conductance 
curves are of very much the same non-linear character, except for the either absence of the humps, or their being 
located at much higher than working-range voltages [0]§. One can see some tendency for developing the humps in 
Fig. 9 (right), and in Fig. 10. Unfortunately, for the stacks we measured, the heating rate was too high to reach 
suggestive voltages presumably corresponding to humps without significant overheating. Perhaps, stacks with linear 
in-plane sizes in the sub-/im range are needed to finally elucidate this issue. 



IV. DISCUSSION 

The previous section in its last part clearly shows that the calculations of the Comment ^[ -> although being 
valid with respect to the general idea of self- heating in the stacks of IJJ, fail to be quantitatively correct. This is 
because the Comment assumed that the resistance is bias independent [ll|, [13 , a simplification which is certainly 
incorrect, see Figs. 9 and 10. It is clear, for instance, that within Newton's cooling scenario, no form of suggested 
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FIG. 10: HgBr2-Bi2212 stacks. The "heating- free" tunnehng-conductance curves deduced out from the pool of numerous "3D" 
I — V — Ts data points taken at about 30 bath temperatures in the range 4 — 90 K. The cr(V) = dl /dV{V) -points in these 
plots correspond to temperatures only different to within ±1.5 K from the temperatures shown along with the curves. Note 
that a{V)(V) oc \V\^, with a > 1 for all temperatures, in contrast to Bi2201-Lacc stacks, for which a < 1, see Fig. 9. Lines are 
guides to the eye. 



bias-independent Rc{T) can account for the observed cr{V) oc < 1, (which is further supported by break- and 

point-contact- junction experiments), [T^. [2(1 [2l| . One has to take into account the highly non-linear /(F)-curves of 
Fig. 9 or 10 to fit better experimental results with help of suggested Newton's cooling scenario. 

Anyway, the self-overheating turns out to be a severe problem in the intrinsic tunneling spectroscopy of supercon- 
ducting gap and pseudogap, mainly due to a very bad thermal conductivity of the host material and large values of 
the gaps involved. 

This problem, although not obviously, can also be acute in some low-resistance break-junction experiments. In- 
deed, the break-junction contact in high-Tc-related experiments is, in overwhelming majority of such experiments, of 
unknown shape, area and location inside the break. Thermal-conductivity properties of such a contact are therefore 
unknown, too. It may happen, that the contact is located at a distance from, say, a metal needle used to produce the 
break 0. Then, the only escape path for the developed Joule heat goes through the poorly conducting Bi2212, i.e. 
one has the same prerequisites for the overheating as in the intrinsic-tunneling experiments. 

Let us make an estimation of the overheating for the (almost) worst case, when the heat is only removed by the ther- 
mally conducting helium gas. We will ignore somewhat better conducting Bi2212 single crystals, leads, gas convection, 
etc. Then, in the spherical geometry, and taking into account the temperature dependence of the thermal conductivity 
of He [2^, which can be approximated by a simple power- low dependence: /^(T) ^ 0.153 (T/300)^*^^^[W/m K], we 
have an approximate estimation: JT[K] ~ 140 • ((5[niW]/r[/im])^-^^, i.e. 140 K for Q = 2 mW and the linear sizes r 
of about 2 /im. 

This formula gives an-order-of-magnitude-correct estimate of the overheating in the case of intrinsic junctions. It 
gives anything between 100 and 1000 K for break-junctions, if we take the range of junction resistances 10 — 3000^7 
and characteristic voltages ~200 mV, and assume that the tunneling junction size is of the order of 100 — 1000 A. The 
latter is most uncertain, but it is very hard to expect that a torn break can have flat and parallel tunneling-j unction 
surfaces of linear size larger than 0.1 /im, unless it is made of fragments of the intrinsic junctions themselves. 



V. CONCLUSIONS 



In conclusion, we have measured directly local overheating in the stacks of intrinsic tunneling junctions exactly 
at places where the bias current is injected. Different Bi-based materials, single-layer, double- layer , an d intercalated 
double-layer ones were investigated. The measurements were motivated by the recent Comment, [Til, [l^ which, by 
using simple Newton's cooling scenario, was able to reproduce high-bias features (pseudogap humps) of our recent 
work ^3 5 thus doubting the validity of those. We confirm the main conclusion of the Comment, that the "pseudogap 
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features" reported in our paper are an artifact of Joule heating. 

However, our measurements have shown that, despite the severe heating in the stacks, there is yet a way of extracting 
the genuine isothermal I-V's. Such /(F, T = const) appear to be highly non-linear in the whole temperature range 
and show notable differences with regard to single- or double-layer compounds. It is also concluded that at moderate 
biases, where the hysteretic multi-branch I-V's are seen, the heating can be very small, on the level of few Kelvin. 

Comparison of intrinsic-junction geometry with other tunneling-spectroscopy techniques, like break-junction ones, 
has revealed, that the heating issues can be acute there as well. All this means that results of tunneling spectroscopy 
on high- Tc materials should be taken with discretion. 
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[0] fit seems to be important to deposit Bi thin film used for calibration and the main thin film in one and the same run. Bi-Ag thermocouple 

deposited in another run apparently had 30% lower thermo-e.m.f. 
[0] JThe number of IJJ is only approximately known to be about 7—10 for that sample. The uncertainty is because the branches were not 

well seen. When the individual "gaps" are not equal (say, due to different doping levels), I-V's show rather chaotic switching between 

different combinations of "good", and "bad" junctions that makes counting them difficult. 
[0] §There are no minima in Rc(T < 300 K) for the samples shown in Fig. 9. This is noted to be inconsistent with their nominal La-content 

of 0.2 (overdoped material, left) and 0.4 (optimal doping, right). The samples may have lost some oxygen during fabrication of the 

stacks. It is therefore not possible to conclude whether th e n on-linearities of the heating- free dl /dV(V)^s shown in Figs. 9 should vanish 

at the positions of minima in the corresponding Rc{Tys |l7j|. 



